Testing single-photon wave packets 
by Hong-Ou-Mandel interference 

Piotr Kolenderski and Konrad Banaszek 



OO 

O 

o 

(N 
C 



43 ■ 
a,; 
^— > ■ 
c ■ 
ctf ■ 
=5 

cr 



> 
in 

(N 

O 

OO 
O 



X 



Abstract — We discuss characterization of single-photon wave 
packets by measuring Hong-Ou-Mandel interference with a weak 
coherent pulse. A complete multimode calculation is presented 
and effects of multiphoton terms in the coherent field as well as 
the impact of source and detection imperfections are discussed. 

Index Terms — Two-photon interference, single photons, quan- 
tum information processing 



I. Introduction 

Single photons are elementary building blocks in optical 
implementations of quantum information processing, commu- 
nication and cryptography [l]-[3]. Apart from the simplest 
protocols that encode and process information using just one 
photon, more elaborate schemes require interactions between 
photons. Experimentally, the most accessible way to realize 
such interactions is to resort to multiphoton interference oc- 
curring in linear optical networks. Despite the limited range 
of transformations that can be achieved this way, supple- 
menting linear optics with auxiliary sources and feed-forward 
operations yields universal quantum information processing 
capabilities [4]. 

Realization of multiphoton interactions by means of linear 
optics places stringent requirements on the modal character- 
istics of single photons. Interference effects require spatio- 
temporal indistinguishability of the single-photon wave pack- 
ets at the input of the linear-optics networks. A sufficient 
test for this condition is the Hong-Ou-Mandel interference 
effect on a 50/50 beam splitter [5], when two incident indistin- 
guishable photons are always found in the same output path. 
One possible realization of such a test is to use two photons 
originating from independent generation events [6]. However, 
more insight into the modal structure of the generated photon 
can be gained by measuring two-photon interference with a 
suitably prepared weak coherent state. If such a state is picked 
off of a stronger coherent beam, its modal characteristics can 
be determined by measuring the remaining macroscopic part 
with one of many classical methods. This indirectly yields 
information about the single-photon wave packet, provided 
that the two-photon interference has been observed. 

The purpose of this paper is to present a full multi-mode 
analysis of Hong-Ou-Mandel interference between a single 
photon and a coherent state. This allows us to address practical 
issues, such as the optimal amplitude of the employed coherent 
state: for its very small amplitude, the counting times would 
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become prohibitively long; on the other hand, increasing the 
amplitude of the coherent state would enhance its two-photon 
and higher terms masking Hong-Ou-Mandel interference. Our 
calculation takes into account also other factors affecting the 
trade-off, such as imperfect preparation efficiency of single 
photons, non-unit detection efficiency and dark counts of the 
detectors. This allows us to link the two-photon interference 
visibility to the actual modal overlap between the single- 
photon wave packet and the coherent state, and determine the 
systematic error of such a measurement. 

Nonclassical interference between a single photon and a 
weak coherent pulse has been observed by Rarity at al. [7], 
who also presented a single-mode model. A measurement 
of two-photon interference with a suitably chosen collection 
of coherent pulse modes enabled the reconstruction of the 
complete single-photon density matrix in the spectral domain 
[8]. Here we extend previously used models to the fully multi- 
mode case and include the complete photon statistics of the 
coherent field. 

II. Beam splitter transformation 

We will consider here the spectral degree of freedom, 
corresponding to optical fields confined to single-mode fi- 
bres. A generalization including spatial degrees of freedom is 
straightforward. We assume that the single-photon wave packet 
is prepared with a probability p, and is described by a density 
matrix g(u),u>'). Thus the quantum state of the field entering 
the beam splitter reads: 

g a = (l-p)|0)<0| +p J dw J duj'g(uj,uj r )^(uj)\0)(0\a{uj') 

(1) 

where a(ui) and a* (to) are the annihilation and the creation 
operators for a frequency u> at the corresponding input port 
of the beam splitter. The coherent state is prepared in a mode 
characterized by a spectral amplitude u(uj), assumed to be 
normalized to one. The associated annihilation operator is 
therefore given by a superposition: 



b = I duj u* (oj)b(uj), 



(2) 



where b(ui) are annihilation operators of monochromatic 
modes entering the second input port of the beam splitter. 
We will use the following notation for a coherent state with 
an amplitude (3 generated in the mode u(u>): 



|/3u(w)>6 = D b (P)\0) = 



exp 



+ / dw/?u(w)&t(w) ) |0) (3) 
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where £>&(/?) — exp(/3ir _ is the displacement operator. 

The total input density matrix is given by the tensor product 
Qin — Qa® \(3u(uj))b(f3u(uj)\. This expression is transformed 
to the output density matrix g out represented in terms of the 
outgoing modes c(ui) and d(u>) leaving the beam splitter by 
substituting monochromatic field operators according to: 

cjuf) - d{ui) ~ c(uj) + d(u>) 

a{uj) = V2 ' m = V2 ^ 
We assume here a flat spectral characteristics of the beam 
splitter over the relevant frequency range. 

III. Coincidence count rate 

The probability of a click on a binary detector monitoring 
one of the outgoing beams is given by the expectation value 
of the following normally ordered operator: 



M c = : I — £ exp 



^—77 J due 1 (uj)c(u) 



(5) 



for the beam c and analogously for the beam d. In the above 
formula, 1 — £ is the probability of a dark count, and 77 is 
the quantum efficiency, assumed to be uniform across the 
spectrum of the detected fields. The coincidence count rate 
Rc is given by the expectation value 



R c =Tr(£ 0Ut : M c <g> M d :). 



(6) 



This quantity can be conveniently expressed in terms of the 
expectation values of a two-parameter operator: 



Z(Vc,Vd) =■ exp ( - / cL>[r] c c ! (lo)c(lo) + r) d $ (w)d(u}]\ 



(7) 
as 

Rc = l- Z[(Z(r), 0)) + (Z(0, T)))} + e(Z( v , 77)). (8) 

In order to evaluate the expectation value (Z(r) c ,T)d)) = 
Tr[g 0Ut Z (j] c , Wd)]> it will be useful to write the output density 
matrix as: 



gout 

P 
2 



where 



(1 - p)g co h+ 
dus I dus' q{w,w')[c < {lo) 



£>coh 



(3u{uS) \ / f3u(uj) 



V2 



V2 



(3u{u>) \ / Pu(lo) 



V2 



V2 



(9) 



(10) 



describes a pair of coherent states in beams c and d prepared 
in a spectral mode described by the mode function u(lo), with 
equal amplitudes /3/\/2. 

The expectation value (Z (rj c , rjd,)) involves then terms of the 
form Tx[Q mb tf(uj)Z{ri c ,rid)\ and Tr[g coh c t (o;)Z(77 c , r} d )c(u)')], 
where the operators for the beam c and be also substituted 
by the field operators for the beam d. The operator products 
appearing in these expressions can be brought to the normally 
ordered form using the formula valid for any bosonic operator 
c [9]: 



Consequently: 

c{uj)Z{r^ c , rjd) = (1 - Vc) Z(r) c , %)c(o;) (12) 

and 

c(u)Z(r) c ,r)d)ct(uj') = 

= (1 - rjc) Z{r] Cl rid) + (1 - Vc) 2 c f (uj')Z(ri c , i] d )c(uj) 

(13) 

With the normally ordered expressions, the expectation value 
Tr [ f3 out Z( rfc, rjd)] can be easily calculated to be equal to: 

Z{VcVd) = 

(l - \{ric + rid)p + \(Vc ~ mfpWT^J e -('^)M 2 /2 

(14) 

where T is the overlap of the single-photon wave packet with 
the mode function of the probe coherent state: 

T = / duj I duj'u*(uj)g(u},u}')u(u'). 



(15) 

Substituting Eq. ( fl4l i into Eq. ((H) we end up with the 
expression for the coincidence rate as a function of the overlap 
T in the following form: 

R C (T) = 1 - £ (2 - V p + \v 2 Pm 2 T^J e-^' 2 + 

+ C 2 (l-vp)e~ vm2 . (16) 

The coincidence count rate measured as a function of the 
delay between the single-photon wave packet and the coherent 
pulse provides information on the spectral characteristics of 
the single-photon wave packet. 

IV. Interference visibility 

The mode overlap T can be recovered from Eq. ( [161 ) by 
comparing the coincidence count rate with the case when the 
modes of the two incident beams are completely distinguish- 
able, corresponding to T = 0. This can be quantified with the 
help of the visibility of the Hong-Ou-Mandel dip, defined as: 



V = 



Rc(0) - Rc(T) 



(17) 



Rc(0) 

After substitution of Eq. ( fToT ) into Eq. ( fT71 ) we obtain that the 
visibility is proportional to the mode overlap T: 



V = c f T 

where the proportionality constant Cf is given by 

iv 2 pW^ ?j2 



2[1 - £(1 - vp)e-'' m2/2 ] (l - £e-'"l' 5 l 2 /2) 



(18) 



(19) 



: exp( 



(— 77c' c) : c = (1 — v)c '■ exp(— rjerc) 



(11) 



and we will refer to it as the correction factor, dependent on 
the parameters of the setup. In the idealized limit of r] = £ = 1 
and a — > the correction factor is equal to one. In the first 
step, it is instructive to consider the case of no dark counts 
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£ = 1 and a small coherent amplitude, which allows one to 
apply an expansion in r\ |/3| 2 . One then obtains: 

c ' wi -G^-i)^i a (2o) 

which shows that the correction factor decreases with the 
increasing coherent state effective amplitude. In the left panel 
of Fig. Q] we plot the correction factor Cf and the coincidence 
count rate Rc for T = as a function of ri\f3\ 2 and r\p 
under the assumption of no dark counts £ = 1, Obviously, 
the coherent pulse amplitude needs to be non-zero to register 
any two-photon event at all, which results in a trade-off for 
its value. 

The behavior of the correction factor c/ is qualitatively 
different in the presence of dark counts characterized by 
£ < 1, as seen in the center and right panels of Fig. Q] 
This is because for very small coherent pulse amplitudes the 
interference dip becomes dominated by dark count events, and 
the visibility carries little information about the mode overlap. 
Consequently, for a fixed product r/p, there is a non-zero 
optimal coherent state amplitude that maximizes the correction 
factor Cf. However, the correction factor never reaches the 
unit value in contrast to the regime free of dark counts. Fig. Q] 
shows that for detectors having 1% and 5% of dark counts the 
highest possible value of the correction factor is close to 0.8 
and 0.7 respectively. 

The statistical uncertainty of the determined overlap T is 
a function of both the correction factor, which attenuates the 
dependence of the measured visibility on T, and the count 
rate of coincidence events. In a simple estimate, the relative 
uncertainty of T can be written as: 

AT 1 

oc (21) 

The expression on the right hand side of the above formula 
can be maximized with respect to r]\P\ 2 for fixed values of £ 
and rjp. Results of this procedure are shown in Fig. [2] which 
provides guidelines for selecting the optimal amplitude of the 
probe coherent pulse. It is seen that the number of detected 
photons in the coherent pulse should be of the order of one. 
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Fig. 2. The dependence of the effective coherent state intensity 7?|/3| 2 on the 
dark count probability £ and the overall single photon production probability 
VP- 
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V. Conclusions 

The method of characterizing single-photon wave packets is 
very simple from the conceptual point of view. However, data 
obtained using an auxiliary coherent pulse to implement Hong- 
Ou-Mandel interference need to be corrected for multiphoton 
events. The optimal amplitude of the pulse has been found for 
a realistic range of parameters. 
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